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After a brief review of the motivations for grand unification, I discuss the main challenges facing 
realistic SUSY GUT model building. Achieving doublet-triplet splitting without fine-tuning is 
chief among them. Symmetry breaking should occur consistently without unwanted Goldston 
bosons, jj. term of order TeV for the MSSM Higgs fileds should emerge naturally, and realistic 
femiion masses with small quark mixing angles and large lepton mixing angles should be gener- 
ated with some predictivity. Significant progress has been made over the years towards achieving 
these goals in the context of supersymmetric 50(10) GUT. A complete 5(9(10) model is pre- 
sented along this line wherein, somewhat surprisingly, the GUT scale threshold corrections to the 
gauge couplings are found to be small. This results in a predictive scenario for proton lifetime. 
An interesting correlation between the d = 6(p ^ e+Tr") and d ~5 (p ^ vK^) decay amplitudes 
is observed. This class of models predicts that both proton decay modes should be observable 
with an improvement in the current sensitivity by about a factor of five to ten. 
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Motivations for unifying the strong, weak and electromagnetic forces are manyfold ||T| ^, ^. 
The experimental observation that electric charges are quantized ([2proton| = |2eiectron| to better 
than 1 part in 10^^) has a natural explanation in grand unified theories (GUT) owing to their non- 
Abelian nature. The miraculous cancelation of chiral anomalies that occurs among each family 
of the SM fermions has a symmetry-based explanation in GUTs. 50(10) GUT, for example, is 
automatically free of such anomalies [||]. GUTs provide a natural understanding of the quantum 
numbers of quarks and leptons. This point is worth emphasizing further. All quarks and leptons of 
a family, including the right-handed neutrino (v'^) needed for generating small neutrino masses via 
the seesaw mechanism, are organized into a 16-dimensional spinor representation of 50(10), as 
shown in Table 1. The gauge symmetry 50(10) contains five independent internal spins, denoted 
as + or — signs (for spin-up and spin-down) in Table 1. Subject to the condition that the number 
of down spins must be even, there are 16 combinations, which form the irreducible 16 dimensional 
spinor of 50(10). The first three spins denote color charges, while the last two are weak charges. 
There are three independent combinations of color spins, identified as the color degrees of freedom 
{r,b,g). Going top down in each column of Table 1, one sees that in addition, there is a fourth color, 
identified as lepton number [|l]|. Thus quarks and leptons are unified under the GUT symmetry. The 
first and the third columns (and similarly the second and the fourth) are left-right conjugates. Thus 
50(10) contains Parity as part of the gauge symmetry. Furthermore, the same 16 multiplet unifies 
quarks with anti-quarks, and leptons with anti-leptons. In fact, 50(10) symmetry is the maximal 
gauge symmetry that is chiral with sixteen particles (members of one family). Hypercharge (and 
thus electric charge) of each fermion follows from the formula Y = ^r(C) — ^L{W), where r(C) 
is the summation of color spins (first three entries) and r(W) is the sum of weak spins (last two 
entries). Thus Y for the e'' field is Y{e') = ^(3) - ^(-2) = 2. Note that Y (and thus Q) must be 
quantized. Such a simple organization of matter is remarkably beautiful and can be argued as a 
strong hint for GUTs. SUSY GUTs have further empirical support from the observed unification 
of gauge couplings at a high energy scale Mx ~ 10^^ GeV. In Fig. 1, left panel, we demonstrate this 
unity of forces in a fully reaUstic 50(10) SUSY GUT [|]. Another remarkable feature of 50(10) 
GUTs is that the small neutrino masses inferred from neutrino oscillation data suggest the scale of 
new physics (v'^ mass scale) to be Myc ~ 10^'* GeV, which is close to Mx- Myc is inferred from the 
effective neutrino mass operator S^y = LLHuHu/Myc (L is the lepton doublet and //„ is the Higgs 
doublet), using my ~ 0.05 eV and (//„) ~ 246 GeV as inputs. In a class of 50(10) models discussed 
further here, Myc ^ /Mpi ~ lO'^ GeV quite naturally [^]. The decay of v'^ can elegantly explain 
the observed baryon asymmetry of the universe via leptogenesis. Finally, as exemplified later, the 
unification of quarks and leptons into GUT multiplets can be quite powerful in realizing predictive 
frameworks for fermion masses, perhaps in association with flavor symmetries. 
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Table 1: Quantum numbers of quarks and leptons. The first 3 signs refer to color charge, and the last three 
to weak charge. To obtain hypercharge, use Y — ^E(C) — ^'L{W). 

While extremely well motivated, constructing fully realistic SUSY GUTs is not so trivial. 
Chief among the challenges is the so-called doublet-triplet (DT) splitting problem. Even in the 
simplest of GUTs, based on SU{5) gauge symmetry the smallest irreducible representation is 
5 dimensional, which contains the SM Higgs doublet. This means that the Higgs doublet will be 
accompanied by a GUT partner, which is a color triplet scalar (5 = 2 + 3 is the relevant math). 
This state must have a mass of order lO'^ GeV, or else it would lead to rapid proton decay. The 
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DT splitting challenge is to naturally make the doublet component of the 5-plet light (of order 
10^ GeV), while maintaining its color triplet partner superheavy. In minimal SUSY SU (5), this 
is done by extreme fine-tuning, of order one part in lO'^. The relevant superpotential is W = 
5//(A24// +M)5h, where the {5h, 5h) contain the MSSM Higgs doublets (//rf,//„), and where 
(24//) = y.diag(l, 1, 1, —3/2, —3/2) breaks the gauge symmetry down to that of the SM. The 
masses of the color triplet and SU {2)l doublet fields are then mj = XV +M, mu = XV — {'i/2)M. 
One chooses XV and M to be both of order 10^^ GeV, but with the condition XV = (3/2) M + 
1^(10^) GeV, so that the doublet remains light. Such a severe fine-tuning raises questions about 
the naturalness of the model. 

In SUSY 50(10) the situation for DT splitting is much better. The adjoint Higgs A(45) in 
50(10) can acquire a vacuum expectation value (VEV) (A) = ia2 (8 Diag (a, a, a, 0, 0). The 
coupling //(10)A(45)//'(10) of two 10-plets would result in heavy color triplets with massless 
SU {2)l doublets - without fine-tuning [^. There are a variety of issues that need to be addressed. 
Symmetry breaking must be complete without unwanted pseudo-Goldstone bosons, the VEV of 
the adjoint should be stable against higher dimensional operators, unification of gauge couplings 
should be maintained, and the theory should be consistent with proton lifetime limits. All these is- 
sues have been successfully addressed recently [fl] by making use a set of small dimensional Higgs 
fields, as shown in Table 2. An anomalous ^ (1 )a symmetry of possible string origin and a Z2 
symmetry are also used, for stabilizing the doublet mass. In Table 2, kis a positive integer, which 
will be taken to be 5. The superpotential consistent with the symmetries is 
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Table 2: '^{1)a andZ2 charges g; and o),- of the superfield i^,. 



W = MaUA'- + ^ (trA^f + ^trA^ + c( ^ZA + ^CC + cis] C' + C'( ^ZA + ^CC + C2s] C 

+ XiHAH' + ( Xh'SZ'^-^+Xh'Z'') ^-^+l2HCC+—AH'CC' . (1) 

This superpotential consistently breaks the 50(10) gauge symmetry down to the SM symmetry 
in the SUSY limit without generating unwanted Goldstone bosons. The first three terms of W 
induce the VEV for A, the fourth and fifth terms guarantee absence of Goldstones [^, and the 
remaining terms achieve doublet-triplet splitting without fine-tuning [^. At the minimum we 
have (C) = (C) = c, (C') = {€')= 0, with c determined by the Fayet-Iliopoulos term of ^(1)^ 
symmetry. The Higgs doublet mass is zero in the SUSY limit to all orders. Once SYSY breaking 
is turned on, the VEV of the A (45) no longer has the zeros, which are modified to be entries of 
order msusY- This in turn induces /i term of order msusY for the MSSM Higgs fields. In Fig. 
I (left panel), the evolution of the three gauge couplings is displayed, which takes into account 
the threshold corrections of the model. Interestingly, the threshold corrections in the model in 
the 10+ 10 sector (of the SU (5) subgroup) cancel between the matter fields and the gauge boson 
fields. Owing to this cancelation, the model becomes very predictive for proton lifetime. We find a 
correlation between the d = 6 gauge boson mediated p ^ e^'nP and the d = 5 Higgsino-mediated 

p vK+ decay amphtudes: Meff ~ lO^'^GeV • (^) ( tH^) ■ controls p e+7r°, 

while Meff controls p — )• vK^. This is plotted in Fig. 1 (right panel) for varying r = M^/Mx (£ is 
a color octet Higgs field). Also plotted are the current experimental limits from these decays. One 
concludes that both modes should be observable with an improved sensitivity of about five to ten. 
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Figure 1: Left panel: Gauge coupling evolution including threshold corrections. Right panel: Correlations 
between M^n and Mx for = 1.5 TeV, = 130 GeV, and a^iMz) = 0.1176. (a): r = 1/200. (b): 
r = 1/250. (c): r — 1/300. The vertical and horizontal dashed lines correspond to the experimentally 
allowed lowest values of Mx and Meff which arise from Umits on F^' (p e^n'^) and r^^{p ^ vK^). 

Realistic and predictive fermion masses and mixings can be obtained within this framework by 
assuming a flavor Q4 symmetry, under which the first two families of 16 form a doublet. The mass 
matrices for up and down quarks, charged leptons and Majorana neutrinos have the form [Q, ^: 

/ £' 0\ / Ki^eS' + ri' \ /^0 0\ 

Mu = -£' a , Mj,^ = -K^^eS' - r]' K'rf.e'^il ^ + ^d,e£ \,MR = m^lio b a\{2) 
\ a \ J \ a + Kd.e'e 1 / \0 a 1/ 

where = 1 and K"e = 3. There are fewer parameters than observables in this setup, which re- 
sults in predictions. A consistent fit for all masses and mixing parameters is obtained with the 
choice a = 0.0508, e = -0.0188 +0.0333/, £ = 0. 106 + 0.0754/, e' = 1 .56 • lO""^, ^' = -0.00474 + 
0.00177/, 1^22 = O.OWe'^-^' at the GUT scale. Along with central values of charged lepton masses, 
we obtain for the quarks, m„(2 GeV) = 3.55 MeV, mc(mc) = 1.15 GeV, md{2 GeV) = 6.45 MeV, 
m.v(2 GeV) = 137.6MeV, mh{mb) = 4.67 GeV. For the CKM mixings we obtain: |V;„| = 0.225, 
iKfa I =0.0414 ,|V„fe I =0.0034 ,|V'f^ I =0.00878 ,17 =0.334 ,p = 0.12, and thereby sin2j8 = 0.663. 
All these are in a good agreement with experiments. For neutrinos, the Dirac mass matrix is ob- 
tained from Mu by replacing e' — )• — 3£'. With 612 ^ 30'^ and 023 — 43" as inputs, we obtain 
m2/m3 ~ 0.13 and 613 ~ 3.6" as predictions. Such a fit is realized by choosing a = 0.0252e^° °'^', 
b = 1.61 • \0-^e-^-^'^^\ and Mo = 1.89- lO^^ GeV. One sees broad, although not precise, agreement 
with data. The model succeeds in obtaining large neutrino mixings along with small quark mixings. 
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